A New Measurement of the Rare Recay K + — > 7r + /i + /i 
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More than 400 K + — > ty + /i + fi~ events were observed in a rare K + decay experiment at the 
AGS. Normalized to the K + — » tt + tt + tt~ decay, the branching ratio is determined to be (9.22 ± 
0.60(stat) ± 0A9(syst)) x 10 -8 . This branching ratio and the fifi mass spectrum is in very good 
agreement with the measurement of the K + — > 7r + e + e~ decay, but deviates significantly from the 
previous measurement. 
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In the Standard Model, the decay K + — > 7r + /i + /i~ 
(-Kw/i/z) proceeds through the same mechanism as K + — > 
7r + e + e~ (AT xee ). It has been recognized for a long time 
that both decays are dominated by long distance contri- 
butions involving one photon exchange Jy-0] and there- 
fore can be described by a vector interaction form factor. 
Consequently, the ratio of the two decay rates depends 
only on the shape of this form factor. Comparisons of the 
Kirfj.^ and K^ee. measurements are of great interest, not 
only within the framework of the Standard Model, but 
also because differences in the form factors measured sep- 
arately in the two decays would indicate new physics, in 
particular those extensions that involve mass dependent 
couplings. The dilepton mass range below 350 MeV/c 2 
is almost inaccessible by other means. A precise mea- 
surement of the decay amplitude is also essential for fu- 
ture investigations of K^^^ in the context of P- and CP- 
violation through measurements of the /i polarization, as 
suggested by many theoretical calculations ||. 

The iYjree mode has been measured several times 
since the mid-70's |J|-||, most recently with high pre- 
cision by this group jg), which observed a sample of 
10,000 events. This recent measurement firmly estab- 
lished that the decay proceeds through a vector in- 
teraction and provided a precise measurement of the 
form factor. The result can be quoted as, Br (-nee) — 
[2.94±0.05 (stat.)±0.13 (syst.)±0M (model)] x 10" 7 and 
5 = 2.14±0.13 (siai.)±0.15 (syst.) where 8 is a parameter 
in the form factor, fv(z) — fo(l + Sz); z = Mu/m,K with 
I = e or /i pof . The slope of the form factor, which is sig- 
nificantly larger than the prediction of meson dominance 
models or leading order chiral perturbation theory, has 
a direct impact on the prediction of the iGr^ branching 



ratio. 

The first observation of the K^ip, decay was reported 
by the E787 collaboration at the AGS in 1997 @ . The 
data sample consists of 13 fully reconstructed three track 
and 221 partially reconstructed two track events, with an 
estimated background of 2.4 and 25 events, respectively. 
The fifi mass spectrum was not analyzed. The branching 
ratio was measured to be [5.0 ± OA(stat) ± 0.7(syst) ± 
0.6(th)] x 10 -8 . The central value used the form factor 
shape measured by 0, corresponding to 8 = 1.31, and 
the last error came from the theoretical uncertainty of the 
form factor shape. This branching ratio measurement 
was found to be 2.2 a below the expectation from the 
then existing K^ ee measurement J?J. The significance 
of the discrepancy becomes greater when this branching 
ratio is compared to the new K^ee measurement, which 
predicts a model independent K^^^ branching ratio of 
(8. 7 ±0.4) x 10 -8 . In view of the important consequences 
of such a discrepancy, an experiment capable of collecting 
a large K^^^ sample with fully reconstructed kinematics 
and improved systematic errors is called for. 

The experiment reported here was performed at the 
Brookhaven Alternating Gradient Synchrotron (AGS) in 
1997, using the E865 detector, which was designed to 
search for the lepton number violating decay K + — > 
Tt + fi + e~ |l^,[l3|]. The data were collected in a period 
of six weeks. The details of the detector design and per- 
formance are discussed elsewhere [J9||l4|] . The components 
relevant to this study are emphasized here. 

With 2.2xl0 12 protons on target, an unseparated 6 
GeV beam of 1.5 x 10 7 K+ , together with 3 x 10 8 ir+ and 
protons, was produced per 1.6 sec AGS pulse. Fig. [I] is a 
schematic diagram of the detector. Decay products were 



measured in a proportional wire chamber based mag- 
netic spectrometer system, together with scintillating ho- 
doscopes and an electromagnetic calorimeter. Muons 
were identified in a range stack consisting of 24 planes of 
proportional tubes situated between iron absorber plates, 
and two sets of hodoscopes (48 pieces each) situated in 
the middle and at the end of the stack. The hodoscopes 
were segmented as vertical strips, separated in the mid- 
dle. 
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FIG. 1. Plan view of the E865 detector. A K^^^ decay is 
superimposed. 

The primary trigger is formed by the A and D ho- 
doscopes and the electromagnetic calorimeter, where 
three charged particle hits are required. Most of the trig- 
ger rate, which is about 70,000/pulse, comes from acci- 
dentals. This raw trigger rate is reduced by a prescale 
factor of 2000. The most frequent K + decay that gives 
three charged particles in the final state is the decay 
K + — > 7t + it + it~ (K T ), and this decay mode is used for 
normalization. For K^^^ decays, the next trigger level 
requires one muon on each side of the detector. Each 
muon is identified as a spatially correlated coincidence 
between the B and C hodoscope hits. The trigger rate 
at this stage is about 2000/pulse, still dominated by ac- 
cidental coincidences. 

In the off-line analysis, events are required to have 
three reconstructed tracks that come from a common ver- 
tex in the decay volume, a reconstructed kaon momen- 
tum consistent with the beam phase space, and a timing 
spread between the tracks consistent with the resolution. 
Fig. H shows the reconstructed kaon mass for K T , in com- 
parison with the Monte Carlo simulation. The mass res- 
olution is a = 2.2 MeV. The Monte Carlo simulation, 
using the GEANT package, takes into account the detec- 
tor geometry, particle decays and interactions with the 
detector, as well as the independently measured efficien- 
cies for each detector component. It reproduces all the 




0.51 
(GeV) 



FIG. 2. Distribution of the reconstructed mm mass for K T 
events. Points represent data, and the histogram is from a 
Monte Carlo simulation. 

For a muon to reach the C hodoscope at the end of 
the muon stack, its momentum must exceed 1.3 GeV. 
Because of the trigger requirement, only events with the 
jx~ on the left side and the n + on the right side of the de- 
tector are accepted. A muon is required to have deposited 
an energy in the calorimeter consistent with a minimum 
ionizing particle, and to have B and C hodoscope hits and 
sufficient muon chamber hits along the projected trajec- 
tory through the muon stack. At this stage, the acci- 
dental background is virtually eliminated, and the 7r/x/i 
candidates are dominated by the background from K T , 
followed by the ir — > fi v decay. The majority of such 
background events have a 7r/i/i mass much smaller than 
Mk- To significantly suppress this background while 
keeping a high efficiency for the signal, a joint likeli- 
hood function is constructed, using the vertex quality, 
K + phase space, and the x 2 0I the tracks. Fig. y is a 
scatter plot of the joint likelihood vs the 7T///Z invariant 
mass. The cluster of events with high likelihood value 
at the kaon mass are the K n ^^ signal, while the events 
below the kaon mass with low likelihood value are back- 
ground from K T . Monte Carlo simulations show that 
it is the pion decays in the spectrometer magnet that 
cause significant mis-measurements of the particle mo- 
menta and spread the ir/ifj, mass upwards to the signal 
region. The cuts based on the likelihood reduce such 
background. Fig. shows the ir/ifi invariant mass distri- 
bution after requiring the joint likelihood function to be 
greater than -13. The observed background events below 
the kaon peak are consistent with the simulation of the 
background from K T . 
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FIG. 3. Scatter plot of the joint likelihood vs Af^ Mfi for 
-Kttmm candidates. The box indicates the cuts for selecting 
K^w events. 

To subtract the background, we model the K nfi n mass 
spectrum with 



n(x) 



— „Pl+P2X+p 3 X 



PA+Pze 



-(x-x Q ) 2 /(2p 2 6 ) 



(1) 



where P\ 1 P2iP^,Pi-,P5-,P& are parameters to be deter- 



mined, 



M„ 



and Xq is the kaon mass. The first 



two terms represent the background, which is subtracted 
in the signal region. 

The same likelihood function is calculated for the nor- 
malization data sample (K T ), as well as Monte Carlo 
events of K T and K viili decays. To verify that the back- 
ground subtraction procedure is reliable, the cut on the 
likelihood function is varied from -18 to -10. Although 
the number of signal events drops from 535 events to 230 
events, and the background fraction improves from 38% 
to 1%, the branching ratio measurement remains con- 
stant within 7%. The cut of -13 is chosen such that the 
statistical error, including the uncertainty in background 
subtraction, is a minimum. There are 430 events in the 
signal region, with 28 background events. The tt/j/jL mass 
resolution is a = 3.3 MeV. 

To examine the decay mechanism, the cos 9 distribu- 
tion (where is the angle between tt + and /i + in the cen- 
ter of mass of the [i + (J-~ pair), and /i + /i~ invariant mass 
distribution are compared to the expectation, as shown 
in Fig ||. The angular distribution is consistent with the 
decay proceeding through a vector interaction, and the 
form factor parameters describing the Af MM spectrum are 
consistent with those of K^ ee . 
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FIG. 4. The nfi/j, invariant mass distribution after requiring 
the joint likelihood to be greater than -13. The shaded area 
is the background according to the fit. 
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FIG. 5. Angular (left) and fi/j, invariant mass (right) distri- 
butions after background subtraction. The points are data, 
the solid histogram is the result of a Monte Carlo calculation 
with 5 — 2.14, and the dashed histogram with constant form 
factor. 

Normalized to the K T branching ratio, the K^^^ 
branching ratio and the form factor slope parameter are 
determined to be 



BR{khh) = (9.22 ± 0.60) x 10" 
o — /.4o_ a95 



(2) 
(3) 



Errors are statistical only, taking into account the uncer- 
tainty in the background fluctuation. 

The similarity between the K T and K^^^ final state 
makes the normalization very reliable. The systematic 
error comes mainly from the uncertainty in the muon 
identification, and the background subtraction proce- 
dure. 




FIG. 6. The \ 2 contours of fits to the Mm distributions 
of the K vlLy . data (solid), and the K wse data(dashed). The 
solid point is the x 2 minimum (Xmin) of the 7r/n/x fit, and the 
open point is that of the 7ree fit. The contours for each fit 
correspond to \ 2 — Xmin + n, n = 1, 4 or 9. Systematic errors 
are included. 

Table | summarizes the systematic errors. There is suf- 
ficient overlap between the primary trigger and the 7r/x/i 
trigger that the ir/i/i trigger efficiency and primary trig- 
ger prescale factor can be precisely determined. To con- 
firm the detector acceptance and efficiencies, many other 
checks were made. The Dalitz plot parameters of the K T 
events were found to be consistent with the published val- 
ues. The muon identification efficiency was checked with 
the K + — > 7r°/i + j/ decay followed by the 7r° — v e + e~j 
decay. 

The K^nn branching ratio divided by the model inde- 
pendent partial branching ratio of K^ee with M ee > 0.15 
GeV is 0.458 ± 0.043. This is to be compared to the ex- 
pectation of 0.432 when 6 = 2.14 is assumed. A more 
complete comparison is shown in Fig. pi where the x 2 
contours are plotted as a function of (/o,<5) for both de- 
cay modes. The consistency is evident. 



Sources 


(TBr/Br 


K T branching ratio 


0.01 


K T prescale factor 


0.01 


7r/i/x trigger efficiency 


0.01 


B, C hodoscope efficiency 


0.02 


EM energy scale 


0.01 


fi Energy loss before C hodoscope 


0.03 


Background subtraction & cuts 


0.03 


Reconstruction Efficiency 


0.015 


Magnetic field map 


< 0.005 


Total 


0.053 



The two normalization samples, namely, K T for K Viljl 
and K + — > 7r + 7r° followed by 7r° — ► e + e~j for K^ ee , are 
compared and found to be consistent within the uncer- 
tainty. 

Our Kjrn^ branching ratio measurement disagrees with 
that of [O. The discrepancy is 3.3 a. 

To summarize, we have measured the K nfifJ , branch- 
ing ratio to be (9.22 ± 0.60(stat) ± 0A9(syst)) x 10~ 8 . 
Both the decay rate and the form factor shape are found 
to be in good agreement with the expectation based on 
the i^Tree measurements, from which we conclude that 
the mechanisms for these two decay modes are consis- 
tent with being the same, as theoretically expected. 
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